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ABSTRACT 

Weak gravitational lensing due to the inhomogeneous matter distribution in 
the universe is an important systematic uncertainty in the use of standard candles 
in cosmology. There are two different descriptions of weak lensing amplification, 
one uses a local smoothness parameter a, the other uses reduced convergence 
7] = 1 + n/l^minl (where k is convergence). The a description involves Dyer- 
Roeder distance Da(oc\z) (5 = 1 corresponds to a smooth universe); it is simple 
and convenient, and has been used by the community to illustrate the effect of 
weak lensing on point sources such as type la supernovae. Wang (1999) has shown 
that the a description can be made realistic by allowing a to be a local variable, 
the local smoothness parameter. The rj description has been used by Wang, Holz, 
& Munshi (2002) to derive a universal probability distribution (UPDF) for weak 
lensing amplification. In this paper, we bridge the two different descriptions 
of weak lensing amplification by relating the reduced convergence 77 and the 
local smoothness parameter a. We give the variance of a in terms of the matter 
power spectrum, thus providing a quantitative guidance to the use of Dyer-Roeder 
distances in illustrating the effect of weak lensing. The by-products of this work 
include a corrected definition of the reduced convergence, and simple and accurate 
analytical expressions for Da(&\z). Our results should be very useful in studying 
the weak lensing of standard candles. 



1. Introduction 

Type la supernovae (SNe la) are our best candidates for cosmological standard candles 
(Phillips 1993; Riess, Press, and Kirshner 1995), and have lead to the startling discovery 
that the expansion of our universe is accelerating (Schmidt et al. 1998; Riess et al. 1998; 
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Perlmutter et al. 1999). This observational fact can best be explained by the existence of 
dark energy in the universe (Garnavich et al. 1998; Wang 2000b). It is important that fur- 
ther and more ambitious SN la observations are carried out (Wang 2000a; Wang & Lovelace 
2001). There are a number of ambitious SN la programs that are active or planned. The 
CFH Legacy Survey, the ESSENCE project at NOAO, and the proposed Large Synoptic Sur- 
vey Telescope (LSST) (http://www.lsst.org/) and Supernova / Acceleration Probe (SNAP) 
(http://snap.lbl.gov), will give us an impressive number of supernovae. The data from such 
observations must be properly analyzed with good understanding of the systematic uncer- 
tainties. 

Weak gravitational lensing due to the inhomogeneous matter distribution in the uni- 
verse is an important systematic uncertainty in the use of standard candles in cosmology 
(Kantowski, Vaughan, and Branch 1995; Frieman 1997; Wambsganss et al. 1997; Holz & 
Wald 1998; Wang 1999; Valageas 2000; Munshi and Jain 2000; Barber 2000; Premadi et al. 
2001). Both the smoothness parameter a (Kantowski, Vaughan, and Branch 1995; Wang 
1999; Sereno et al. 2001, 2002) and the reduced convergence rj (Valageas 2000; Munshi and 
Jain 2000; Munshi and Wang 2003) have been used in the literature to study the weak lensing 
of standard candles. 

Magnification of a standard candle can be expressed in terms of the ratio of our actual 
distance to the standard candle and our distance to it in a completely smooth universe. Such 
distances, Djy(a\z), are solutions to the Dyer-Roeder equation (Dyer and Roeder 1973; Kan- 
towski 1998), with the matter inhomogeneity parametrized by a local smoothness parameter 
a (Wang 1999). Due to their simplicity and convenience, Dyer-Roeder distances with fixed 
values of a are used to illustrate the effect of weak lensing of point sources, although they do 
not correctly describe weak lensing (Hamana, Martel, and Futamase 2000). However, since 
the Dyer-Roeder equation follows from the general equation for light deflection in general 
relativity (for a given mass density) (Schneider et al. 1992), it is valid when shear can be 
neglected (as is the case for weak lensing amplification), as long as mass density is not fixed 
as in the usual applications of the Dyer-Roeder equation. In Wang (1999), a new and valid 
way of using the Dyer-Roeder equation has been developed: by allowing a dispersion in the 
possible values of a (with probability density function derived from the matter power spec- 
trum) due to the inhomogeneous distribution of matter, we can very efficiently compute the 
lensing effect and intuitively understand the weak lensing of cosmological standard candles. 

The reduced convergence 77 has been used by Wang, Holz, & Munshi (2002) to derive 
a universal probability distribution (UPDF) for weak lensing amplification. To achieve a 
unified understanding of weak lensing magnification, it is important that we understand 
how the local smoothness parameter a and the reduced convergence 77 are related. 
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This paper is organized as follows. We derive analytical solutions to the Dyer-Roeder 
equation in Sec. II, to be used in the rest of the paper. In Sec. Ill, we relate the local 
smoothness parameter a and the reduced convergence 77. In Sec. IV, we derive the variance 
of a in terms of matter power spectrum. Sec.V contains a summary and discussions. The 
Appendix contains simple fitting formulae that comprise an accurate expression of the Dyer- 
Roeder distance. 



2. Angular diameter distance as function of the smoothness parameter a 

In a Hubble diagram of standard candles, one must use distance-redshift relations to 
make theoretical interpretations. The distance-redshift relations depend on the distribution 
of matter in the universe. Wang (1999) has shown that we can define the local or direction 
dependent smoothness parameter a via the Dyer-Roeder equation (Dyer and Roeder 1973); 
there is a unique mapping between a and the magnification of a source. In this section, we 
express the angular diameter distance to a standard candle in terms of a. This will be used 
in the next section to relate a to the reduced convergence 77. 

The local smoothness parameter a essentially represents the amount of matter that 
can cause magnification of a given source. Since matter distribution in our universe is 
inhomogeneous, we can think of our universe as a mosaic of cones centered on the observer, 
each with a different value of a. 

In a smooth Friedmann-Robertson- Walker (FRW) universe, a — 1 in all beams; the 
metric is given by ds 2 = dt 2 — a 2 (t)[dr 2 /(l — kr 2 ) + r 2 (dQ 2 + sin 2 6 d<p 2 )], where a(t) is the 
cosmic scale factor, and k is the global curvature parameter (Q k = 1 — Q m — Q x = ~k/ H^). 
The comoving distance r is given by (Weinberg 1972) 



where "sinn" is defined as sinh if fi^ > 0, and sin if 0^ < 0. If fi^ = 0, the sinn and Sl^'s 
disappear from Eq.(l), leaving only the integral. Qxfx{ z ) is the contribution from dark 
energy; the dimensionless dark energy density f(z) = px(z) / px(z = 0). For a cosmological 
constant, Qx = ^a, and fx{z) = 1- The angular diameter distance is given by (1a(z) = 
r(z)/(l + z), and the luminosity distance is given by d^z) = (1 + z) 2 dA(z). 

However, our universe is clumpy rather than smooth. According to the focusing theorem 
in gravitational lens theory, if there is any shear or matter along a beam connecting a 
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source to an observer, the angular diameter distance of the source from the observer is 
smaller than that which would occur if the source were seen through an empty, shear- 
free cone, provided the affine parameter distance (defined such that its element equals the 
proper distance element at the observer) is the same and the beam has not gone through 
a caustic. An increase of shear or matter density along the beam decreases the angular 
diameter distance and consequently increases the observable flux for given z. (Schneider 
et al. 1992) For studies of weak lensing magnification (with convergence \k\ < 0.2 (Barber 
2000)), we can ignore shear and consider convergence only, which corresponds to the matter 
in the beam. 



For given redshift z, if a mass-fraction a of the matter in the universe is smoothly 
distributed, the largest possible distance for light bundles which have not passed through 
a caustic is given by the Dyer-Roeder distance, the solution to the Dyer-Roeder equation 
(Dyer and Roeder 1973; Kantowski 1998): 



f \ d 
9{z) dz 



9(z) 



dD A 

dz 



+ -&n m (i + z ) 5 D A = o, 



D A (z = 0) = 0, 



dD f 
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2 = 
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g(z) = (l + z) 2 E{z). 



(2) 



The angular diameter distance for a given smoothness parameter a and redshift z, 
Da(6[\z), can be obtained via the numerical integration of Eq.(2). However, it is instructive 
to express D A (a\z) in terms of a. 

Integrating Eq.(2) gives 

g(z) - cH ' 1 = -^a f D A {a\z>) dz', (3) 

dz 2 J g(z') 

where we have used the boundary condition on dD A /dz at z — 0. Integrating the above 
equation gives 

n <~\ \ tj-i r dz> 3 n ~ r dz> r'( i +^) 5 n t ~\ «\ a « <*\ 

D A (a\z) = cH I — - -n m aj o — ^ L-^J- D A (a\z ) dz . (4) 

Now we exchange the order of integration in the double integral above, i.e., J* dz' J* dz" = 
dz" f*„ dz' . Hence 

D ^ ^ c ^ [ w) - t a " a i' w DAan dz " L m (5) 
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Now we have 

D A (&\z) = ~Q m & f" dz' (±±4r- D aW) [X(Z) ~-i Z ' )] + D A (a = 0\z), (6) 
2 Jo g{z') cH 

where the affine parameter 

Hz) = cffo" 1 f 44- (7) 



/o gW 

Note that D A (a = 0\z) = X(z). Wang (2000b) has shown that D A (a = l\z) = r(z)/(l + z), 
where r(z) is the comoving distance in a smooth universe. 

No approximations have been made in obtaining Eq. (6) from the Dyer-Roeder equation 
Eq.(2). We can solve the Eq.(6) perturbatively by replacing D A (a\z') in the integral on the 
right hand side with its previous order approximation. In the first order perturbation, we 
replace D A (a\z') in the integral by its 0th order approximation, D^\a\z') = D A (a = l\z'). 
This gives us 

D ( j\a\z) = D A (a = l\z) aR min (z)+D A (a = 0\z) = D A (a = l\z) [1 + K min (z)(a - 1)] , (8) 

where we have defined 

„ D A (a = 0\z) 3 Q m f z M+zT r(z') 

using D A (a = l\z) = r(z)/(l + z) (Wang 2000b). Note that k min (z) < 0. 
After four more iterations, we find the fifth order perturbative solution 

r, D ( A ^ Z l , = l-\hmin{z)\{a-l) {!-& [C^z) - C 2 (z) a + C,{z)a 2 - C 4 (z)a 3 + 0(a 4 )} } , 
D A {a — l\z) 

(10) 

The functions Ci(z) (i = 1, 2, 3, 4) are given by 

c, (z) fi^A+j^BA^M | S „ (z0 | Ci _ iM [AM-y-)] 

' 2 \k ml „(z)\ J E(z') D A (a= l\z) 1 m ' nK " ' cH^ 1 1 ' 

where we have defined C = 1. 

Fig.l shows the difference between the angular diameter distance given by the solution of 
the Dyer-Roeder equation and our analytical expansion to 4th order in a (see Eq.(10)), with 
^min(z), and Ci(z) (i = 1,2,3,4) computed numerically using Eq.(ll). Clearly, the accuracy 
of Eq.(10) is practically unlimited, as additional terms can be added in a straightforward 
manner as needed. 

Note that we have not yet made any assumptions about dark energy; the results of this 
subsection are valid for all cosmological models. Simple and accurate fitting formulae for 
^min(z), and Ci(z) (i = 1,2,3,4) are given in the Appendix. 
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3. Reduced convergence and the local smoothness parameter a 

Wang, Holz, and Munshi (2002) have derived a universal distribution for the weak 
lensing magnification distribution of standard candles, valid for all cosmological models, 
with arbitrary matter distributions, over all redshifts. Their results are based on a universal 
probability distribution function (UPDF), P{rj), for the reduced convergence, 77. For a 
given cosmological model, the magnification probability distribution, -P(/i), at redshift z is 
related to the UPDF by P(/i) = P(rj)/2 \K min \, where 77 = 1 + (/i - 1) / (2\K min \) , and n min 
(the minimum convergence) can be directly computed from the cosmological parameters. 
The UPDF can be well approximated by a three-parameter stretched Gaussian distribution, 
where the values of the three parameters depend only on £ v , the variance of rj. In short, all 
possible weak lensing probability distributions can be well approximated by a one-parameter 
family. Each alternative cosmological model is then described by a single function £ v (z)- This 
method gives P(fJ-) in excellent agreement with numerical ray-tracing and three-dimensional 
shear matrix calculations, and provide numerical fits for three representative models (SCDM, 
ACDM, and OCDM). 

It is of interest to relate the local smoothness parameter a to the reduced convergence 
77, as both can, and have been used to parametrize weak lensing of standard candles. 

The magnification factor of a source is given precisely by (Schneider et al. 1992) 

(l — K,y — y 

where k is convergence (which depends only on matter within the beam), and 7 is shear 
(which depends on matter outside of the beam), k and 7 can be computed from the deflection 
potential. Ignoring shear, we can relate magnification /i to convergence k: 

" = (nbp (13) 

which is valid for \k\ < 0.2 (equivalent to /i < 1.56) for all cosmological models at all redshifts 
(Barber 2000). 

In terms of angular diameter distances, magnification can be defined as (Schneider et al. 
1992) 



D A (a = 1) 2 



D A (a) 

Comparing Eqs.(13) and (14), we have 



(14) 



k=1-j^L- y (15) 
D A (a = 1) 
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The minimum convergence corresponds to the minimum magnification p m i n , which occurs 
when a = (empty beam, when the smoothly distributed matter fraction in the beam is 
zero), i.e., p = 0, or 5 = (p — ~p)/~p — — 1. Since the minimum magnification can only be 
achieved when 7 = (see Eq.(12)), the true minimum convergence is given precisely by 

ft ■ 1 ^(* = °) (16) 
D A (a = 1) 

as we have defined in the previous section. Note that the smallest possible magnification, 
Umin, corresponds to k min by definition, 



Pmin 



D A (a = l) 2 



D A (a = 0) 



>,1 K"min) 



(17) 



In the context of weak lensing, k is the projected density contrast, and n m in is the 
projected density contrast when 5 — — 1 along the line of sight from observer to source 
(Valageas 2000; Munshi and Jain 2000): 



K r . 



Thus 



2(cH 1 yJ r( X ) 
dX = rW = Jgsinn (VWx) ■ (18) 

ft^ = tidz'%fr(z')(l + z) [X(z)-X(z')] 

where we have used Eq.(9). Fig.2 shows the difference between K min {z) and K min (z) for three 
representative cosmological models. 

To see why n min differs from k min , let us write them in more intuitive forms: 

Kmin(z) = --——j dz — — — - D A (a = 1 z , z) , (20) 

2 cif Jo ^(^ ) D A (a = l\z) 

where we have used D A (a = l\z) = r(z)/(l + z) (Wang 2000b), and D A (a = 0\z',z) = 
(1 + z') [X(z) — X(z')} is the angular diameter distance of a source at z from a fictitious 
observer at z' for a = (Schneider et al. 1992). Similarly, we can write 

mm( j " 2^J dz "eRT ^(5 = i|.) TT^ ' (21) 

where r (x(z) — x(- 2 ')) /(I + z ) * s the angular diameter distance between z and z' in a com- 
pletely smooth universe (i.e., filled beam, (5 = 1). Since the true minimum convergence 
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should correspond to the minimum magnification, which occurs when the beam connecting 
the source and observer is empty (i.e., a — 0), k min (z) is the true minimum convergence, 
while K m i n {z) is its approximation when the filled beam (instead of the empty beam) angular 
diameter distance is used along the line of sight from the observer to the source. 

K m in(z) has been mistakenly associated with /i min in the literature. Note that \K m i n \ < 
\K m in\, hence 1/(1 — n m in) 2 underestimates the minimum magnification fi m i n . Some of the 
consequences of the subtle difference between k min and K min will be presented in Sec. IV. 

Fig. 2 shows that the correction to K m i n (z) required for the accurate calculation of the 
minimum convergence K m i n (z) is quite modest (less than 2% at z < 1) in cosmological models 
that fit current observational data (the ACDM model). 

Let us define a modified reduced convergence given by 

I K-min | 

which should be compared with the reduced convergence defined by Valageas (2000) and 
Munshi and Jain (2000), 

V = 1 + r^T (23) 



i] is related to rj by 



Since 



rj = l + [fj-l]p™!\. 1 (24) 



D A {a = 0) - D A {&) 

D A (a = 1) ' (25) 
we find that 

fj = a {1 - (a - 1) [d(z) - C 2 {z) a + C 3 {z)& 2 - C A (z)a z + C(a 4 )] } = 77(a), (26) 
where we have used Eq.(10). 

In the limit of weak lensing, both the modified reduced convergence 77 and the reduced 
convergence r\ approach the local smoothness parameter a: r\ ~ fj ~ a. 



4. Dispersion in the local smoothness parameter due to inhomogeneous 

matter distribution 

It is most convenient to illustrate the effect of gravitational lensing by comparing the 
angular diameter distance in a smooth univerve (which is equal to Dyer-Roeder distance 
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with smoothness parameter a = 1 (Wang 1999)) with Dyer-Roeder distances with a ^ 1. 
Since this is often done by the community, it is important that we qualify such illustration 
by noting the following: 

(1) A realistic description of weak lensing of point sources using the smoothness parameter 
a requires generalizing a to be a local parameter, the local smoothness parameter (Wang 
1999). 

(2) The distribution of the possible values of the local smoothness parameter a is directly 
related to matter distribution in our universe. 

(3) The variance of a depends on the matter power spectrum. It provides a quantitative 
guidance in the use of Dyer-Roeder distances to illustrate the effect of gravitational lensing 
of point sources. For example, instead of arbitrary values of a, one should use a values which 
are one or two standard deviations away from the mean value of (a) = 1. 

Next, we derive the variance of a in terms of the matter power spectrum. This is done 
through the relation of a to the reduced convergence rj derived in Sec. III. 

In Sec. Ill, we have shown that K min is the true minimum convergence and corresponds 
to the minimum magnification p m i n , while K min is the approximation to k min in the weak 
lesning limit. This suggests that we need to use fj = 1 + n/\k min \, instead of rj = 1 + K,/\K min \, 
as the reduced convergence. Subsequently, we modify the formulae in Valageas (2000) as 
follows. 

The convergence k (Bernardeau, Waerbeke, and Mellier 1997; Kaiser 1998) is modified 
to be 

«=|n m cff - 1 £dz^^-S, (27) 
where z s is the source redshift, 5 = (p — p)/p, and 

w(z, z s ) = 2 ^ (1 + z f{l + z s ) [X(z s ) - X(z)} . (28) 

The minimum of k occurs when p = 0. It is straightforward to check that the minimum 
of the k as given in Eq.(27) is h min (see Eq.(9)), which does correspond to the minimum 
magnification (see Eq.(17)). 

The variance of fj is similar to that of rj (Valageas 2000), and given by 

^-{—) cH ° I iz W) fewi" (29) 

with k min (z) given by Eq.(9), and 

W^n^f ^^W*(Vk9 ), (30) 



where A 2 (k,z) = 4irk 3 P(k, z), k is the wavenumber, and P{k,z) is the matter power spec- 
trum. The window function W(Vk0 ) = 2J 1 (Vk9 ) / (Vk9 ) for smoothing angle 6 . Here J\ 
is the Bessel function of order 1. 



Since 77 and a are related (see Eq.(26)), we have 



1 N 

1 



da 



dfj 



(31) 



where is the variance of a. Hence 




drj 



(32) 



9a 



with 77(a) given by Eq.(26), and ^ given by Eq.(29). 



5. Summary and discussions 



We have related the local smoothness parameter a to the reduced convergence r] = 
1 + K/\K min \. This establishes the connection between two different approached to the weak 
lensing of point sources. 

For simplicity and convenience, it is likely that the community (especially observers) will 
continue to use Dyer-Roeder distances for various values of the smoothness parameter a to 
illustrate the effect of weak lensing on point sources. In order to provide meaningful guidance 
to such illustrations, we have devived the variance of the local smoothness parameter a in 
terms of the matter power spectrum (see Eqs.(32), (29), and (26)). 

Wang, Holz, & Munshi (2002) have used rj to derive a universal probability distribution 
(UPDF) for weak lensing amplification. We propose a corrected definition for the reduced 
convergence: 77 — 1 + K/\k min \, with fj, min = 1/(1 — k min ) 2 . We have shown that the true 
minimum convergence is R m in (defined as the convergence corresponding to the minimum 
magnification fj, m i n ), while K m i n (z) (as commonly defined in weak lensing literature) is only its 
approximation when the filled beam (instead of the empty beam) angular diameter distance 
is used along the line of sight from the observer to the source. 

Note that n m i n does not correspond to /i m j„. However, the correction needed for n m in 
is modest; K m i n is smaller than k m in by less than 2% in a fl m = 0.3 universe at z < 1 (see 
Fig. 2). Thus we don't expect dramatic corrections to the various published weak lensing 
statistics. 
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As a by-product of our derivation relating the local smoothness parameter a and the re- 
duced convergence fj, we obtained simple and accurate analytical expressions for the angular 
diameter distance, D A (a\z), which are solutions to the Dyer-Roeder equation. We find an 
expression of D A (a\z) as a polynomial in the local smoothness parameter a (with coefficients 
given by integrals) that can be made arbitrarily accurate, and is valid for all cosmological 
models (see Eq.(lO) in Sec. 2). In the Appendix, we give simple analytical approximations 
of the coefficients in Eq.(lO) for flat models with a cosmological constant, and open models 
without dark energy. These analytical approximations give angular diameter distances that 
have an accuracy better than 0.1% up to a redshift of 2. 

Our results should be useful in studying the weak lensing of standard candles. 

We thank Ron Kantowski for helpful discussions, This work was supported in part 
by NSF CAREER grant AST-0094335, and the NSF REU program at the University of 
Oklahoma. 



A. Analytical Approximations for K min (z) and Ci(z) 

Eq.(10) gives the Dyer-Roeder distance in terms of a polynomial in the local smoothness 
parameter a. For clarity and convenience, let us now derive analytical approximations for 
K m i n (z), and Cj(z) (i = 1, 2, 3, 4), for the special case of the dark energy being a cosmological 
constant. 

First, let us expand the integral expressions in Eqs.(9) and (11) in the limit of small z. 
We find for z <C 1, 

Kmin(z) - ~^fz 2 [1-0.5 (1 + 2q ) z] , 

d(z) ~ ^mZ 2 [l-0.5(l + 22o)2], 

d{z) ~ Cln 2 m z 2i [l-0.5i (l + 2q )z], (* = 2,3,4) (Al) 
where C 2 = j^j, C 3 = j^o, C 4 = 394 2 400 , and q = Vt m /2 — Q A . 
Ci(z) (i = 1, 2, 3, 4) can be rewritten as 

Ci(z) = C w (z) \k min (z)\ , 

d(z) = CMdWd-^z), (^ = 2,3,4) (A2) 
where Cio (i = 1, 2, 3, 4) are nearly constant. We find that 

C w (z)~0.3, C 20 (z)~^, C 30 (z)~^-, C 40 (z)^^-. (A3) 
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Guided by the small z limit, we are able to find simple analytical forms to approximate 
ii m in(z) and Aq(z). 

For a flat universe, we found: 

Q m z 2 

Krnin{z) = ^ [di In (c^ 3 + 1) + 1] , 

ai(fi m ) = .06354 fi- L5498 + .26 88^ 8586 + 1.58 12, 

a 2 (ft m ) = 3.6049 (^ 18 V 1082 ^ - fi m ) - .1751, 

a 3 (On) = n^ 306 e - 6022n - - 1.4427 ft m + . 5472. (A4) 

and 

C 10 (z) = 0.3 + (-0.00548 - 0.02725) In [E(z) - z (1 + g )] • (A5) 

For an open universe, we find 

O ^ 2 (1 + 0.5*) f 0.19 \n[E(z)-z(l + q )] \ 

K ™^ Z >- 4 \ X (^1.08 Z 0.45(|lnn m |+0.04r 7 /' ^ 

and 

C w {z) = 0.3 + [fi^ 52 In (ft™ 557 - 0.284) + 0.302] In [E{z] - z (1 + q )\ . (A7) 

Figs. 3-4 show the accuracy of our analytical approximations for K m i n (z) (see Eqs.(A4)(A6)) 
and Cio(z) (see Eqs.(A5)(A7)) as function of fl m , for redshifts up to 2 and 5 respectively. 
Each figure displays the maximum of the ratio of the difference between our analytical ap- 
proximation and the exact expression over the latter. 

Fig.5(a) shows the difference between the angular diameter distance given by the nu- 
merical solution of the Dyer-Roeder equation and our analytical expansion to 4th order in a 
(see Eq.(10)), with K min (z) and Ci(z) (i = 1,2,3,4) given by Eqs.(A3)-(A7). The accuracy 
shown in Fig. 5 (a) is primarily limited by our approximations for k min (z), Eqs.(A4) and (A6). 

Fig.5(b) is the same as Fig.5(a), except with k m i n (z) given by Eq.(9), with the empty 
and full beam distances, D A (a — Q\z) and D A (a = l|z), computed by integrating the Dyer- 
Roeder equation. 

Comparison of Figs. 3-4 with Fig. 5 shows that the errors in estimated distances using 
analytical approximations of Eqs.(A3)-(A7) are much smaller than the errors in K m i n (z) and 
A (z). This is because 



AD A (a\z) A\k min (z)\ 
D A (a\z) 



D A (a = l\z) _ 



D A (a\z) 



^ A\K min (z)\ ^ Ag ^ 
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The last inequality arises from the fact that D A (ct — l\z) and Da{ol\z) differ by less than 
~ 20% for a < 2. 

For all practical purposes in the weak lensing of standard candles, the local smoothness 
parameter a < 2, since the probability for a > 2 is vanishingly small for all cosmo logical 
models at all redshifts (Wang 1999). Note that the probability distribution of a becomes 
narrower and peaked closer to a = 1 as redshift z increases, as the universe is smoother 
on the average at high z (Wang 1999). Our analytical approximations, given by Eqs.(10), 
(A3)-(A7), give an accuracy of better than about 0.1% in all cosmological for a < 2 (see 
Fig.5(a)). If the exact expression for k m i n (z) (see Eq.(9)), instead of Eqs.(A4) and (A6), is 
used, the accuracy in the estimated distance is better than 0.02% for all cosmological models 
for a < 2 (see Fig. 5(b)). 
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Fig. 1. — The difference between the angular diameter distance given by the solution of the 
Dyer-Roeder equation and our analytical expansion to 4th order in a (see Eq.(lO)). 
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Fig. 2. — The difference between n min {z) and n min (z) for three representative cosmological 
models. 
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Fig. 3. — The accuracy of our analytical approximations for K min (z) (see Eqs.(A4)(A6)). 
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Fig. 4. — The accuracy of our analytical approximations for Ci (z) (see Eqs.(A5)(A7)). 
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Fig. 5. — (a) The difference between the angular diameter distance given by the solution 
of the Dyer-Roeder equation and our analytical expansion to 4th order in a (see Eq.(lO)), 
with K m i n (z), and Ci(z) (i = 1,2,3,4) given by Eqs.(A3)-(A7). (b) The same as (a), except 
with K m i n (z) given by Eq.(9), with the empty and full beam distances, Da{o. = 0,z) and 
Da(oi = l,z), computed by integrating the Dyer-Roeder equation. 
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